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Abstract: This study analyzed nucleotide sequences from the mitochondrial cytochrome oxidase submit (COI) gene 
region (450 bp) to investigate the genetic structure of the oriental river prawn ( Macrobrachium nipponense ) among nine 
populations from the Yangtze and Lancang Rivers. A total of 79 individuals were collected for this work. Eighty-nine nu- 
cleotides were found to be variable, resulting in 46 haplotypes. Among the nine populations, the population from Kunming 
shows the greatest level of variability (hk = 1.000) x =0.028), whereas the population from Chongqing exhibits the lowest 
level of variability (hk = 0.700, x =0.008). Analysis of molecular variance suggested that of the total genetic diversity» 
9.66% was attributable to inter-population diversity and the remainder (90.34%) to differences within populations. A 
molecular phylogenetic tree constructed using the Neighbor-joining (NJ) method showed that the 46 haplotypes were as- 
signed to two clades associated with geographic regions. These results provide basic information for the conservation and 
sustainable exploitation of this species . 
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The oriental river prawn © Macrobrachium nippo- tion naturally in most regions of China. These facts» 
nense ) is one of the few Macrobrachium species that combined with its popularity with Chinese consumers» 
has a temperate distribution and is capable of reproduc- make it an attractive aquaculture species ( New» 
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2006). Over the past 20 years M . nipponense has be- 
come a commercially important aquaculture species. 
However» M . nipponense farming has recently substan- 
tially declined due to disease and the small size of the 
cultured animals (Fu et al, 2004a). The wild popula- 
tion has undergone a steady decline at the same time» 
mainly because of overexploitation, pollution, and con- 
struction of hydropower projects. However, the genetic 
variation among M .nipponense populations is largely 
unknown, and there is little information on population 
differentiation. Information on genetic variation is vital 
for the design and implementation of adequate manage- 
ment strategies for the species. Furthermore, in order 
to relieve pressure on wild stock, and to improve in- 
dustry products through selective breeding, develop- 
ment of a domestication program for M .nipponense is 
important, requiring baseline information concerning 
genetic background information and genetic variations . 

Until recently» all major genetic studies on this 
prawn have been based on morphological characteris- 
tics» cytogenetics» interspecific hybridization and poly- 
ploidy breeding (Mashiko» 1983, 1990, 1992; Mashiko 
& Numachi> 1993; Qiu et al» 1994, 1997; Mashiko & 
Numachib, 2000; Fu et al» 2004a; b; Zhao» 2006a)+ 
with the exception of Ni ©2003) and Zhao et al 
(2006b), who simply studied genetic structure of this 
prawn. There has been no research on population varia- 
tion in addition to the studies mentioned above. There- 
fore» it is necessary to study the genetic diversity of this 
species» especially using DNA tools. 

As a genetic marker» mitchondrial DNA (mtD- 
NAD has a few advantages over other markers. This in- 
cludes the virtue of its maternal» nonrecombining mode 
of inheritance, its rapid pace of evolutions and exten- 
sive intraspecific polymorphism (Avise et al» 1987). 
mtDNA is widely used in the study of genetic variation 
in organisms, including some crustaceans. The object 
of this study is to investigate the genetic structure of the 


oriental river prawn CM . nipponense > within nine popu- 
lations from southwest China» using nucleotide se- 
quences from the mitochondrial cytochrome oxidase 
submit I (COI) gene region» and to provide baseline 
information for the conservation and sustainable utiliza- 
tion of the wild genetic resources of the oriental river 
prawn. 


1 Material and Methods 


1.1 Sample and DNA extraction 

A total of 79 individuals of Macrobrachium nippo- 
nense were collected from nine localities. Details of 
each locality are shown in Tab 1. Samples were pre- 
served in 75% ethanol until DNA was isolated for ge- 
netic analyses. 

Ventral muscle tissue was dissected and incubated 
in a standard buffer [0.06 mmol/L ethylenedi- 
aminetetraacetic acid (EDTA), 0.1 mmol/L Tris (pH 
8.6) 0.5% sodium dodecyl sulphate (sps>] overnight 
at 37°C in the presence of Proteinase-K. After buffer 
incubations DNA was extracted using standard phenol/ 
chloroform extraction protocols and concentrated using 
column purification. The purified DNA was dried» dis- 
solved in TE buffer and stored at -20C. 

1.2 DNA amplification 

Primers were designed using Primer Premier 5.00 
(Singh et al» 1998) based on the complete mtDNA of 
M . rosenbergii Miller et als 2005). Segments of the 
COI gene for each individual were amplified using PCR 
of COI F (5'- TIT ATC TTC GGA GCG TGA GC -3°) 
and COI R (5'- AGT TAT TCC TGG GGC TCG TAT G 
-3'). The PCR reaction mixture of 30 pl containing 1.2 
units of Tag DNA Polymerase» 1 x reaction buffers 5 
mmol/L MgCl, 0.4 pmol/L primer, 250 pmol/L of 
each dNTP; and up to SOng of genomic DNA. The PCR 
cycling included pre-denaturing for 3 min at 94°C and 
30 cycles of 1 min at 94°C 1 min at 50°C, and 1 min 
at 72°C; followed by final extension for 5 min at 72°C. 


Tab. 1 Population codes, locations, sample size, and genetic diversity of the nine populations 





Haplotype diversity Nucleotide diversity 


Population 


Locality 


Sample size 


Latitude 


Longitude 





h) (x) 

sZ Suzhou» Jiangsu 9 31°10'N 120°38'E 0.917 0.013 

NC Nanchang, Jiangsu 9 29°22'N 116°20'E 0.806 0.016 

HZ Huzhou, Zhejiang 10 30°40'N 120°29'E 0.889 0.016 

cQ Chongqing 5 29°32'N 105°45'E 0.700 0.008 

YD Yidu, Hubei 10 30°38'N 111°45'E 0.800 0.016 

WH Wuhan» Hubei 9 30°33'N 114°19'E 1.000 0.016 
SM Simao» Yunnan 10 22°40'N 101°27'E 0.933 0.017 

KM Kunming» Yunan 25%00'N 102°40'E 1.000 0.028 

XSBN Xishuangbanna, Yunan 21°10'N 100°40°E 0.861 0.027 
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Samples wese sent to Shanghai DNA BioTchnologies 
Co» Ltd. for sequencing. For each sample, sequenc- 
ing was performed in both directions. 
1.3 Data analysis 

Sequence chromatograms were viewed and edited 
manually using Chromas 2.13. Once edited; multiple 
alignments were performed using Clustal X 1.81 
(Thompson et al» 1997). Haplotype diversity Ch) 
(Nei & Tajima, 1981) and nucleotide diversity Cz > 
(Nei, 1987) were calculated using Dnasp 4.0 (Rozas 
et aly 2003). An analysis of molecular variance 
CAMOVA; Excoffier et al» 1992) was performed to 
partition the total phenotypic variance into intra- and 
inter-population variance using ARLEQUIN 3.01 CEx- 
coffier et al» 2005). The neighbor-joining © NJ > 
method was used to reconstruct the phylogentic rela- 
tionships among haplotypes with MEGA (version 2.1: 
Kumar et al» 2001). 


2 Results 


2.1 Sequence variation and haplotype distribution 

A total of 450 bp mtDNA COI sequences were am- 
plified successfully from 79 individuals of nine popula- 
tions» resulting in the identification of 46 unique haplo- 
types defined by 89 variable sites (GenBank accession 
No: EF102435 - 102480). The mean total nucleotide 
composition was A = 30.7%» C= 21.2%» G = 19.2%; 
T = 28.9%. Among the 46 haplotypes» seven (15.2%) 
were shared by different populations; other haplotypes 
were unique to the corresponding population. Four 
haplotypes ChapS, hap8-» hapl1 and hap15) were 
shared between two of the populations» and three other 
haplotypes Chap4, hap10 and hap21) were shared by 
another three of the populations (Tab. 2). 
2.2 Population genetic diversity and genetic di- 

vergence 

The amount of COI gene sequence variation among 
the nine populations is summarized in Tab. 1. The Kun- 
ming (KM) population exhibits the greatest level of vari- 
ability Ch = 1.000, æ = 0.028), whereas the Chongqing 
CCO? population exhibits the lowest level of variability Ch 
= 0.800, x = 0.008). AMOVA identified that a high 
proportion of the total genetic variance was attributable to 
variations among populations (90.34%), whereas only 
9.66% occurred within populations (Tab. 3). 

The fixation index € Fşr) and genetic distance 
Nei’ s measure) between pairs of populations are 
shown in Tab. 4. These data suggested that the 


Huzhou (HZ> and Suzhou (SZ) populations were ge- 
netically most similar © Fsr = - 0.055), while the 
Yidu (YD?) and Chongging CCQ? populations were 
most dissimilar  Fşr = 0.241). 
2.3 Phylogeny reconstruction 

The neighbour-joining tree (Fig. 1) was constructed 
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Fig. 1 Neighbor-joining tree of COI haplotypes 


of Macrobrahcium nipponense 
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Tab. 2 Distribution of the 46 haplotypes in Macrobrchium nipponense populations 











Haplotype 

Population 1 2 3 4 

1234567890123456789012345678901234567890123456 
Suzhou 3; l uf 11 
Nanchang 4 2 1 1 1 
Huzhou 23 2 1 E 
Chongqing at 1 3 
Yidu 1 41 3 Ñ 
Wuhan ai 11 1 Lote tet. i 
‘Xishuangbanna 2 2 1 1 
Simao 1 1 1 1 1 1 al 
Kunming Ttt 1 1 1 af 





with the complete data set of 46 haplotypes. Two main 
clades were identified. Samples from the Wuhan 
(WH), Yidu CYD), Nanchang (NC); Huzhou (HZ> 


Tab. 3 Analysis of molecular variance (AMOVA) for 
the Macrobrachium nipponense populations 











and Suzhou (SZ) populations formed one cluster mid- Source of Degree of Variance Percentage of 
‘ es Z ny 
dle and lower reaches of Yangtze River). The other i aaa - =— r ee 
$ A mong populations i r 
clade was composed of the Chongging (CQ), Kunming Within populations 70 0.443 90.34 
(KM?, Simao SM) and Xishuangbanna (XSBN? pop- Total 78 0.490 100 
ulations» which were collected from upriver reaches of Fer: 0.097 
Yangtze River and Lancang River. 
Tab. 4 Fixation index ( Fsr) (below) and genetic distance (above) in nine 
populations of Macrobrachium nipponense 
Kunming Suzhou Nanchang Huzhou Chongqing Yidu Wuhan Xishuangbanna Simao 
Kunming ae 0.041 0.033 0.040 0.019 0.030 0.031 0.031 0.025 
Suzhou 0.042 = 0.020 0.015 0.029 0.024 0.021 0.031 0.038 
Nanchang 0.099 0.094 = 0.019 0.020 0.018 0.017 0.027 0.030 
Huzhou 0.057 -0.055 0.069 = 0.028 0.023 0.020 0.031 0.037 
Chongqing 0.134 0.177 0.204 0.190 = 0.018 0.018 0.021 0.015 
Yidu 0. 104 0.143 0.197 0.156 0.241 = 0.020 0.027 0.026 
Wuhan 0.000 0.030 0.097 0.035 0.131 0.102 = 0.027 0.029 
‘Xishuangbanna 0.030 0.111 0.167 0.125 0.208 0.151 0.069 e 0.025 
Simao =0,003 0.075 0.120 0.089 0.165 0.097 0.034 0.003 





3 Discussion 


This study used nucleotide sequences from the 
COI gene region to investigate the genetic structure of 
Macrobrachium nipponense among nine populations 
with a total of 79 individuals. The percentage of A + T 
base composition (59.6%) was much higher than C + 
G» which coincides with other crustacean protein-cod- 
ing genes (Miller et al» 2005). 

The present study of sequencing an approximately 
450 bp fragment of this gene region revealed 46 haplo- 
types based on polymorphisms at 89 nucleotide sites . 
Cook et al (2002) investigated the genetic structure of 
throughout the 
Queensland in Australia. Their study revealed that se- 


M . australiense rivers of western 


quences of a 607 bp fragment of the COI gene revealed 
17 haplotypes in 28 individuals from six sites» based 
on polymorphisms at 36 nucleotide sites. Bruyn et al 
(2004) studied the phylogeography of M . rosenbergii 
from the Lake Carpentaris region. They amplified the 
602 bp COI gene for all samples taken from 378 indi- 
viduals, resulting in 43 unique haplotypes defined by 
59 variable sites. Compared with the above studies» 
the present research resulted in more variable sites and 
haplotypes due to the larger collection range and more 
informative characters of the COI gene. 

Pairwise Fsy analysis and Nei’ s genetic distance 
analysis indicated that the different of 
M . nipponense between the Huzhou (HZ?) and Suzhou 


genetic 


(SZ) populations is the smallest. which corresponds 
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with the geographic distance. In this study» the Kun- 
ming (KM? population presents the greatest level of 
variability Ch = 1.000, x = 0.028), while the 
Chongqing ©CQ> population exhibits the lowest level of 
variability Ch = 0.800, x = 0.008). The smaller (5 
individuals) sample size might be part of the reason for 
the lower variability observed in the Chongqing CCQ?) 
population. Furthermore, the development of small hy- 
dropower stations in this region have had adverse eco- 
logical effects on the aquatic system of the dammed 
streams, and have inevitably affected lower reaches, 
causing habitat fragmentation within dammed rivers and 
significant habitat changes downstream (Zhou, 2004). 
This would result in a decrease of effective population 
size and a decline in genetic diversity in the Chongqing 
CCQ) population. 

The AMOVA analysis provided corroborating evi- 
dence for the genetic structure obtained from Nei’s ge- 
netic diversity statistics. Of the total genetic diversity, 
9.66% was attributable to inter-population diversity 
and the remainder 90.34%) to differences within 
populations. This indicates that there is a small differ- 
entiation among the nine populations, due to the high 
dispersal ability of this prawn and sample locality of 
this study. The Fsr value (0.097) suggested that al- 
though the differentiation between the nine populations 
is small» there is moderate genetic divergence among 
populations . 

The mtDNA data resolved two divergent mono- 
phyletic clades» with the first representing the middle 


and lower reaches of Yangtze Rivers and the second 
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